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Abstract. This work focuses on the properties of dusty tori in active galactic 
nuclei (AGN) derived from the comparison of SDSS type 1 quasars with mid- 
Infrared (MIR) counterparts and a new, detailed torus model. The infrared data 
were taken by the Spitzer Wide-area InfraRed Extragalactic (SWIRE) Survey. 
Basic model parameters are constraint, such as the density law of the graphite 
and silicate grains, the torus size and its opening angle. A whole variety of 
optical depths is supported. The favoured models are those with decreasing 
density with distance from the centre, while there is no clear tendency as to the 
covering factor, ie small, medium and large covering factors are almost equally 
distributed. Based on the models that better describe the observed SEDs, prop- 
erties such as the accretion luminosity, the mass of dust, the inner to outer radius 
ratio and the hydrogen column density are computed. 



1. Torus Model 

The torus geometry adopted to describe the shape and the spatial distribution 
of dust is the so-called flared disk, that is a sphere with the polar cones removed. 
Rs size is defined by its outer radius, Rout, and the opening angle, of the torus 
itself. The dust components that dominate both the absorption and the emis- 
sion of radiation are graphite and silicate. The location of the inner radius, 
Rin, depends both on the sublimation temperature of the dust grains (1500 and 
1000 K, for graphite and silicate, respectively) and on the strength of the accre- 
tion lu minosit y. We adopted the absorption and scattering coefficients given by 
iLaor Drain (|l993l ^ for dust grains of differen t dimensions, weighted with the 
standard MRN distribution ( Mathis et al-llToT^ . Grains dimensions range from 



0.005 to 0.25 micron for graphite, and 0.025 to 0.25 micron for silicate. The gas 
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density within the torus is modeled in such a way to allow a gradient along both 
the radial and the angular coordinates. 

The central source is assumed to be point-like and its emission isotropic. Its 
spectral energy distribution is defined by means of a composition of power laws 
with different values for the spectra index in the UV, optical and IR. A numerical 
approach - the A-iteration method - was adopted to solve the radiative transfer 
equation. A geometrical grid is defined along the three spatial coordinates, and 
the main physical quantities (dust density and temperatures, electromagnetic 
emission, optical depth, etc.) are computed with respect to the center of the 
volume elements defined by the grid. When computing the total incoming energy 
on a given volume element a very accurate computation of the optical depth, is 
done. Furthermore, since silicate grains have a lower sublimation temperature 
with respect to graphite, we account for the fact that the innermost regions of 
the torus contain only graphite and are silicate free 

Finally, the global SED is computed at different angles of the line-of-sight 
with respect to the torus equatorial plane, in order to account for both type 1 
and type 2 objects emission and includes three contributions: emission from the 
AGN, thermal emission and scattering emission by dust in eac h volume element . 
For an example of an emitted spectrum see Fig. ^ We refer to lFritz et al.l ( 20061 ) 
for a detailed description of the model developed to compute the emission of dust 
in AGN. 




A[/x] 

Figure 1. Emission spectra as a function of wavelength for 10 different 
lines-of-sight inclinations from = 0° (edge-on; lower curve) to 'J = 90° 
(face-on, upper curve) at regular steps of 10°, for a geometrical configuration 
with r(9.7) = 10.0, Rout/Rin = 300, torus opening angle 9 = 140° and a 
density that depends on both the angle and the distance from the centre, like 
p{r,e) oc r-ie-6|™"WI 
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2. Data 

The MIR data used here are taken from the SWIRE ELAIS Nl and N2 and the 
Lockman fields, and were obtained from February 2004 through July 2004, with 
both IRAC and MIPS. For the purposes of this work only the four IRAC bands 
and MIPS 24 micron are used. The SWIRE catalogues we use we re processed by 
the SWIRE collaboration. Details about the data can be found in lLonsdale et al.l 
( 2004^ .ISurace et al.1 (j200.^ and Shupe et al., m prep. 

The Sloan Digital Sky Survey has validated and made publicly available 
its Data Release 4 (DR4), covering the entire Lockman and SWIRE EN2 fields 
and a part of the SWIRE ENl field. A total of 280 spectroscopically con- 
firmed quasars lie within the SWIRE fields covered by the SDSS DR4 spectro- 
scopic release. Their i-band magnitudes reach 19.1 for objects with redshifts 
typically less than 2 .3 and go up to a magnitude deeper for higher redshifts 
( Richards et"Il^l2002^. For a detailed analysis of the properties of the sample in 
ENl see lHatziminaoglou et al.l (|2005l ^ 

3. First Results 
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Figure 2. Examples of SED fits for three AGN of our sample. The fluxes 
are given in Jy, the wavelengths in A. 

Each observed SED is compared to a total of 720 models, with 10 lines 
of sight each (from an equatorial to a pole-on view). Examples of best fits are 
shown in Fig. [21 Here we summarize some of the results of the fitting procedure: 

• Size of Torus: Overall, Rout/R-in and the torus covering factor are equally 
distributed between the allowed values (30, 100, and 300; 95%, 75% and 
50%, respectively). However, both inner and out radii clearly increase with 
z, whereas larger covering factors models tend to fit realtively low redshift 
{z < 1) objects. 

• Density: For ~70% of the objects the density (written as p (r, 9) oc r^e~^^^°^^^^^) 
decreases radially from the centre while for another ~25% it remains con- 
stant. 

• Optical depth: This work addresses, among other things, the possibility of 
low optical depth tori. The SED fitting resulted in some ~45% of objects 
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Figure 3. Left panel: Outer torus radius, Rout, as a function of redshift, 
z, showing a tendency of larger tori at higher z. Right panel: Mass of dust, 
MDust, as a function of z, indicating more massive tori at higher redshifts. 



with Tg.y > 1.0 and '^55% with rg.y < 1.0. This impHes that, for some 
configurations, type 1 AGN could be seen even when the line of sight 
intercepts the torus. 

The distribution of outer radius of the torus, Rout; and the mass of dust, MDust, 
with z is shown in the left and right panels of Fig. |31 respectively. Moust is 
computed summing the individual sample elements of the best-fit model. The 
increase of both Rout and M^ust with redshift could indicate larger and more 
massive tori at higher redshifts, however it could also simply be the manifestation 
of the Malmquist bias, as both quantities scale with the accretion luminosity, 
Lqcci an input of our model that is proportional to the bolometric luminosity 
computed from the data. Note that Moust does not refer to the total mass of 
the torus, that can be obtained by adding the mass of gas, typically '^100 times 
larger than that of the dust. 
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